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Abstract 

We consider and compare four Hamiltonian formulations of thermostated mechanics, three of 
them kinetic, and the other one configurational. Though all four approaches "work" at equilib- 
rium, their application to many-body nonequilibrium simulations can fail to provide a proper flow 
of heat. All the Hamiltonian formulations considered here are applied to the same prototypical 
two-temperature "0 4 " model of a heat-conducting chain. This model incorporates nearest-neighbor 
Hooke's-Law interactions plus a quartic tethering potential. Physically correct results, obtained 
with the isokinetic Gaussian and Nose-Hoover thermostats, are compared with two other Hamil- 
tonian results. The latter results, based on constrained Hamiltonian thermostats, fail correctly to 
model the flow of heat. 

PACS numbers: 05.70.Ln, 05.45.-a, 05.45.Df, 02.70.Ns 
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I. INTRODUCTION 



Computational "Thermostats" arose as a means for controlling numerical simulations 
of both equilibrium and nonequilibrium stationary states. Without thermostats systems 
driven away from equilibrium heat up. With thermostats the heat generated by irreversible 
processes can be steadily removed, making it possible to simulate nonequilibrium steady 
states Because most Hamiltonian-based mechanics problems conserve energy, novel 

nonHamiltonian ideas are typically required when thermostats are to be included. Nev- 
ertheless, several approaches to Hamiltonian thermostats have been developed. Here we 
consider various approaches pioneered by Bill Ashurst, Carl Dettmann, Denis Evans, Wm. 



G. Hoover, 
25 years [ 



bm Leete, Gary Morriss, Shuichi Nose, and Les Woodcock over a span of about 



er, lorn Leete, Uary IV. 
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The simplest thermostat type maintains a (nearly) constant kinetic energy by "rescaling" 
the velocities at the end of each computational timestep. For a more elegant, but equivalent, 
continuous approach to rescaling, see the "Gaussian isokinetic" method described in Section 
lib and illustrated in Section VII. For # Cartesian degrees of freedom the kinetic-theory 
relation, 



= * mg2 = *jf_ kT\ 



kin 



2 ^ 2m 2 ' 

defines the kinetic temperature Tki n . Thermostats based on can be applied to an entire 

system, or separately to two or more subsystems. With this ad hoc rescaling approach, 



p^ p yJ(K/K ) , 

there is no difficulty in maintaining several temperatures in specified parts of a single 
nonequilibrium system, as did Ashurst in his Ph D thesis research 11| . 

Kinetic theory and the Gibbs-Boltzmann development of thermodynamics based on ideal- 
gas thermometry certainly suggest that the kinetic temperature, T^ Q , is both the simplest 
and the most fundamental of the many possible "temperature" types. There are alternatives. 
Among the many possible temperatures based on particle coordinates, rather than velocities, 
the simplest — on which we focus here — is based on the mean-squared forces, (F 2 ). This 
definition of a configurational temperature is derived directly from the potential energy $ 
and its space derivatives: 

kT = <(W) 2 )/(V 2 $> = -(F 2 )/(VF> . 



This somewhat cumbersome relation appears as an aside in Landau and Lifshitz' classic 



1958 text 
discussed 



12] . It was also recently rediscovered by Rugh, and has subsequently been much 
131 ]. This configurational representation of temperature follows most simply from 



an integration by parts, in the canonical ensemble: 

J dqV 2 § = V$ ; {-d/dq)e- q,/kT = (V$/A;T)e-* /fcT — > 

J dqV 2 $e-*' kT = (1/kT) J dq(V$) 2 e-^ kT . 

About ten years after Woodcock and Ashurst's introduction of the rescaling isokinetic 
thermostat in the early 1970s, Hoover and Evans discovered that the isokinetic velocity- 
rescaling equations of motion can be derived from Gauss' Principle of Least Constraint 5] . 
Fifteen years after this discovery, Dettmann and Morrissj^] found a straightforward Hamilto- 
nian approach to these same motion equations. Previously, Leete and Hoover had derived a 
Hamiltonian formulation which likewise maintained the velocity-based kinetic energy K({q}) 
constant [a 

A Hamiltonian generating the Gauss'-Principle isokinetic equations of motion was not 
found until 1996 [4]. In the meantime Nose had discovered his canonical-ensemble ther- 
mostated dynamics, which is related to two rather different Hamiltoniansj?], ^]. The purpose 
of the present work is twofold. First, we explore the relation of the Hoover-Leete kinetic 
thermostat work to its successors. A 2007 literature search shows only a single reference to 
We also explore the consequences of a "Landau-Lifshitz" configurational thermostat 
based on the canonical-ensemble definition rediscovered by Rugh. 



II. THE NOSE-HOOVER THERMOSTAT 

The Nose-Hoover equations of motion [l|, 0] are a standard approach to carrying out 
canonical-ensemble dynamics at a particular temperature To. This approach uses feedback, 
controlling the fluctuating kinetic energy with a "friction coefficient" or "control variable" 
(. The frictional feedback force is — (p. The time rate of change of the friction coefficient 
is based on the current value of the kinetic energy K(t) relative to the desired mean value 
Kq(T ). The equations of motion also include a characteristic time r which determines the 
time-range over which the feedback acts: 

{q = p/ m; p = F - Cp}; { = [(K/K ) - l]/r 2 . 
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In most cases, including our 4 model studies described in detail in Sections V-IX, a rea- 
sonable choice of r corresponds to a typical collision time or vibration time. Such a choice 
typically provides Gibbs' equilibrium canonical phase-space distribution with a Gaussian 
distribution for the friction coefficient £: 

f(q, P ,() = e-*/ kT e- K / kT e-#^ 2 / 2 

Evidently the amplitude of the fluctuations of ( varies as ^l/^r 2 , where # is the number 
of degrees of freedom. Thus its effect on the dynamics vanishes in the large-system limit. 

Nonequilibrium situations are quite different. By using two or more different friction 
coefficients (or "thermostat variables") temperature differences can be established, leading 
to heat flow. Then the corresponding phase-space distributions are no longer smooth and 
Gibbsian, but instead become fractal, with the underlying phase-space trajectories satisfying 
the Second Law of Thermodynamics infinitely more probable than those violating the LawQ]. 
In the equilibrium case there is a close connection of Nose- Hoover mechanics to Hamiltonian 
mechanics. We describe this connection next. 



A. The Dettmann-Nose Hamiltonian 



In 1984 Shuichi Nose discovered a Hamiltonian consistent with Gibbs' constant- 
temperature canonical ensemble, 

f{q,p) oc exp(-H Nos6 /kT) . 

He was able to derive the Nose-Hoover equations of motion described above from his Hamil- 
tonian by an artificial "time scaling" . 

In July of 1996 Carl DettmannQ discovered a simpler approach to the Nose-Hoover 
equations of motion. Dettmann's vanishing Hamiltonian is 

n2 



T^Dcttmann = s ^Nose = / , ~ h S 

2ms 



. 



Provided that Dettmann's Hamiltonian is set equal to the special value, zero, Nose's time- 
scaling variable s can be eliminated. Then the Nose-Hoover equations of motion result 
without the need for any time scaling. It needs to be emphasized that this TC — trick 
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does not work for the prototypical nonequilibrium case of a system with two different tem- 
peratures. In a two-temperature system the two different variables SHot/Coid are necessarily 
coupled and cannot both be eliminated, so that Dettmann's 7i = trick does not work. 

There is a special case of Nose- Hoover dynamics that deserves special mention, the "isoki- 
netic" case in which the temperature is constant. This corresponds to the choice r — > 0. 
For isokinetic mechanics Dettmann and Morris discovered the special Hamiltonian detailed 
next. 



B. The Dettmann-Morriss Isokinetic Hamiltonian 

Soon after Dettmann discovered a Hamiltonian route to the Nose-Hoover motion equa- 

n 

tions, Dettmann and Morriss[:4|] discovered a related approach to isokinetic (constant tem- 
perature, provided that temperature is defined by the kinetic energy) Hamiltonian dynamics: 

H DM = Ke + ^ 2K ° - K e-^ 2Ko = ; K = V ^- . 

2m 

As usual, the equations of motion follow by differentiation: 

{ q = dH/dp = pe+* /2K ° ; p = -dH/dq = Fe^ 12 ^ } . 
The accelerations which follow, using the identity, 

K/K = e-* /Ko , 

are then exactly the same as those from the isokinetic equations of motion given above: 
{ q = (F/m) -p^iF- q/2K Q )e + ^ 2 ^ = (F/m) - (q } ; C = £ F ■ q/(2K ) . 

III. THE HOOVER-LEETE ISOKINETIC THERMOSTAT IN ITS LAGRANGIAN 
AND HAMILTONIAN FORMS 

The Hamiltonians corresponding to the isokinetic Gauss'-Principle and Nose-Hoover 
velocity-based approaches are not the only such means of thermostating equilibrium sys- 
tems. A straightforward application of nonholonomic Lagrangian mechanics, as well as the 
familiar Hamiltonian mechanics outlined in Leete's Master's thesis jfj, lead to another type of 
isokinetic mechanics. These two equivalent forms of mechanics, both of which can be used to 



keep the velocity-based kinetic energy K(q) constant, while allowing the momentum-based 
kinetic energy K{p) to vary, proceed by modifying the {q}, rather than the {p}, equations. 
To distinguish the two different kinetic energies we use the notation: 

K{4) = -£?f;K<p) = ■££.. 

Let us begin with the Lagrangian case: 

£hl = K(q) - $ + \[K(q) - K ] , 

where K{q) is the velocity-based kinetic energy, J2 y<? 2 , and the Lagrange multiplier A has 
the task of maintaining K(q) at its initial value, K , as the motion proceeds. Lagrange's 
equations of motion follow from the usual textbook differentiations of the Lagrangian with 
respect to the velocities and coordinates: 

{ P = ~~^T7~~ = mq(l + A) ; p — mq(l + A) + mq\ = *^ HL = F } . 
oq oq 

Now multiply the p equation by q and sum: 

mq(l + A)g + m ^ = + 2K o^ = J2F-q = -^ . 
The value of the Lagrange multiplier A follows: 



($o - $) = y/4K(p)K - K = 2K (1 + A) - 2K 
2K ~ 2K ~ 2K ' 

where the last expression, which completes the identity, follows from the Hamiltonian given 

just below. 

With Leete's help, Hoover discovered, in 1979, that the velocity-based kinetic energy K(q) 
can alternatively be kept constant by using the (constant in time) Hoover- Leete Hamiltonian: 

Hhl = - £hl = \/4K K(p) + $ -K = H = % + K , 

where $ is again the usual potential energy and $ is its initial value. 

The constancy of the kinetic energy is easy to see. The equations of motion, 



OH . OH p / K . 



imply that the velocity-based kinetic energy does not vary: 



K . 
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At the same time, except in the equilibrium case with a large number of degrees of freedom, 
there is no guarantee that the momentum-based kinetic energy K(p) is similar in magnitude 
to K . But, provided that { mq = p } — > K(p) = K(q) = K initially, Leete's equations of 
motion do match the Gaussian isokinetic ones to second order, 

r •■ F EF-q 

{q = — (i \ ; C = ^ , 

through the second derivatives, { q }, but differ in the third derivatives, where the mean- 
squared constraint force in the isokinetic case is less than the Hoover-Leete analog, in ac- 

nn 

cordance with Gauss' Principle of Least Constraint LSI] . 

The constraint of constant temperature in Hamiltonian mechanics contradicts the ther- 
modynamic notion that the energy and temperature cannot be varied independently in a 
system of fixed composition and volume. Both E and T are constants of the motion, us- 
ing Leete's approach. Provided that the initial conditions are wisely chosen, with E and 
T corresponding to the same thermodynamic state, this approach can certainly be used to 
determine equilibrium properties. We will see, in Section VIII, what the (rather strange) 
consequences of this thermostat are away from equilibrium. 



IV. LANDAU-LIFSHITZ' CONFIGURATION AL THERMOSTAT 

A (much) more complicated Hamiltonian thermostat, conserving the force-based configu- 
rational temperature T con , can be based on straightforward (though quite tedious) holonomic 
Hamiltonian mechanics. A proper configurationally thermostated Hamiltonian conserves not 
only the Hamiltonian but also the configurational temperature, 

N N 

Here k is Boltzmann's constant, which we set equal to unity in the numerical work. A 
molecular dynamics simulation based on this definition of temperature follows standard 
Hamiltonian mechanics, as the temperature constraint is just a (complicated) holonomic 
(coordinates only) constraint. The simplest procedure begins with the system Lagrangian, 
augmented with a Lagrange multiplier A which constrains { T con , T con , T con }. The corre- 
sponding Hamiltonian has the form, 

H = H{q}) + K({p}) + A(T con - T ) . 
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For success, the initial conditions have to be chosen to correspond to both the desired 
temperature, T con = T , and to the condition T con = 0. Then, to begin the analytic work, 
differentiate the temperature equation, T con ({g}) = T , twice with respect to time. The two 
time differentiations, using the chain rule, give first {q} and then {q}. By substituting the 
constrained equations of motion: 

{ q\ = Si = pi = Fi({q}) - AViT con } , 

for the {q}, we obtain the Lagrange multiplier A. For further details of this calculation, see 
Section IX. The numerical work can be checked by noting that both T = T con and 7i are 
constants of the motion when the calculation is error free. 



V. AOKI AND KUSNEZOV'S <^ 4 MODEL FOR HEAT CONDUCTION 

The (f) 4 model gets its name from the functional form of a quartic "tethering potential" , 
^Teth; which links each particle to a fixed lattice site with a cubic restoring force: 

$Teth = \^~J ~* FTcth = ' 

In the pedagogical simulations which follow we will choose the tethering force constant equal 
to unity, k — 1. Nearest-neighbor particle pair interactions in the <p A model are governed by 
Hooke's Law: 

$Pa ir = £W2)(M-rf) 2 - 
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Aoki and Kusnezov have carried out a variety of pedagogical heatflow simulations for this 
model in both one and three space dimensions [fj], [l5]. The 4 model has a finite Fourier 
conductivity in one dimension. The numerical work carried out by Aoki and Kusnezov 
established the temperature dependence of the thermal conductivity, 3/T 4 / 3 . 

For simplicity in the one-dimensional work which we carry out here, we always choose the 
nearest-neighbor separation d equal to unity. We also choose the strength of the Hooke's-Law 
interaction fi equal to unity. As was abundantly demonstrated by Aoki and Kusnezov, the 
combination of a site-based tethering potential with a Hooke's-Law pair potential provides 
the usual Fourier conductivity, with the heat flux proportional to the (sufficiently small) 
temperature gradient. In what follows we apply four different thermostat constraints to 



S 



this potential model, using four different dynamical approaches: Nose-Hoover, Gaussian 
isokinetic, Hoover-Leete isokinetic, and Landau-Lifshitz isoconfigurational. 

A numerical solution of the heat flow equation for the one-dimensional chain, 

f = V(3T~ 4/3 VT) ± aT , 

can be obtained by Runge-Kutta integration, with a rescaling of the temperatures within 
the reservoir regions at the end of each timestep accounting for the source and sink terms 
±aT in the flow equation. Figure 1 shows a fully-converged temperature profile obtained 
in this way, using hot and cold temperatures of 0.26 and 0.24, with 400 mesh points. The 
spatial gradient operations were replaced by finite-difference approximations: 

(VT), = (T i+1 - T^O/2 , 

with appropriate subscript changes accounting for the periodic boundary conditions. 



VI. RESULTS WITH THE NOSE-HOOVER THERMOSTAT 

With a periodic system composed of four parts, "Hot", "Newtoni", "Cold", and 
"Newton 2 ", each part containing N/4 particles, the equations of motion are 

{ mq = F — ( Hot mq } ; (uot = [(K/K Hot ) - l]/r 2 ; 

{ mq = F — (coidmq } ; Ccoid = [(K/K Co id) ~ 1]/V 2 5 
in the hot and cold regions, and 

{ mq — F } 

in the two Newtonian regions. The 4 force Fj depends upon the three coordinates 
{xi-i, Xi, Xi + i} and also includes the tethering force, -Fr e th = —k{xi — x ) 3 , where xq is the 
lattice site for the ith particle in the perfect lattice. The initial velocities and displacements 
were selected randomly within ranges: 

— A„ < {v} < +A V ; -A, < {5q} < +A q , 

with A^ or A q chosen to reproduce the desired stationary reservoir temperatures, 0.26 and 
0.24. The average temperatures, both kinetic and configurational, from the last half of a 
400-particle, 40-million timestep simulation, 

< t < 100, 000 (equilibration) ; 100,000 < t < 200, 000 (averaging) , 

9 



are shown in Figure 2, along with a comparison solution of the macroscopic heat flow equa- 
tion. The agreement is quite good, as would be expected from Aoki and Kusnezov's work. 

VII. RESULTS WITH THE GAUSSIAN ISOKINETIC THERMOSTAT 

Simulations with two Gaussian isokinetic thermostats, with the same timestep, dt = 
0.005, and the same equilibration and averaging intervals, produced similar results, main- 
taining the initial kinetic temperatures throughout. The averaged profiles, along with a 
corresponding continuum solution, are shown in Figure 3. Again the agreement, for both 
the kinetic and the configurational temperatures, is perfectly satisfactory. 

The instantaneous (at time t) values of the total heat added (from time to time t) to 
the "hot" reservoir as well as that removed from the "cold" reservoir are shown in Figure 
4. The sum of the two necessarily oscillates about zero as t increases. Straight lines drawn 
through the data provide an estimate for the Fourier conductivity in agreement with the 
results of Aoki and Kusnezov's work. The two conventional thermostat types, Nose- Hoover 
and Gaussian isokinetic, show that either method can be used to simulate the simple two- 
temperature heat flow problem. 

VIII. RESULTS WITH THE HOOVER-LEETE HAMILTON! AN ISOKINETIC 
THERMOSTAT 

Results with the Hamiltonian isokinetic thermostats: 

r . P I -^Hot/Cold . TP/ \ ^ 

{9Hot / coid= mV~^r ;p= 

where the fixed hot and cold kinetic temperatures (mq 2 /2/c)Hot/Coid are fixed by the ini- 
tial conditions, but come to differ from the unconstrained momentum-based temperatures 
{p 2 /2mk) Hot /Cold- Unlike the Nose-Hoover and Gaussian isokinetic profiles, the Hamiltonian- 
based profiles show no temperature gradients at all. See Figure 5. Instead there are dis- 
continuities between the fixed kinetic temperatures and the adjacent constant values of 
the unconstrained Newtonian regions. In the latter regions the kinetic and configurational 
temperatures equilibrate, and match the unconstrained configurational temperature of the 
thermostated regions. In the thermostated regions there is no such equilibration. 
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Evidently there are no heat fluxes in the system either. By constraining not only the 
energy (through the constant Hamiltonian driving the motion) but also the kinetic tem- 
peratures, the system is evidently prevented from supporting heat flow. The Hamiltonian 
nature of the motion equations also prevents the formation of the multifractal phase-space 
distributions associated with nonequilibrium stationary states l|, |2j. A more detailed analy- 
sis of the reasons for the failure of this Hamiltonian approach cou 
detailed analysis of the microscopic equations describing heat flux 



d perhaps be based on a 



IK- 



IX. RESULTS WITH THE LANDAU-LIFSHITZ HAMILTONIAN ISOCONFIGU- 
RATIONAL THERMOSTAT 

For the special 4 model considered here, with d = 1, the force Fi is 

Fi = -k8* + /i[5 i+1 - 25i + <Ji_i] . 
The pair-potential contribution to V 2 7i is 2/x, for the two nearest-neighbor interactions, 

d 2 



■(li/2)(\xi-x i±1 \-l) = n . 



The tethering potential contribution is 3nSf = V 2 (nSf/A): 



N N 

The time derivative of the configurational temperature is a quotient of products of single- 
particle sums: 



£(2/i + 3/^ 2 ) [E(2/i + 3K5 2 )] 2 ' 
where 

F = /i[v i+1 - 2vi + Vi-t] - 3n5 2 Vi . 
Another time differentiation gives the second derivative: 

2J2FF + 2J2F 2 72 J2F 2 [J2Kv8f 



-^Newton , — , / „ t-o \ " ' 



E(2/x + 3^ 2 ) [£(2/i + 3/^ 2 )] 3 

24E^E^ + 6Ei r2 (^E[a^ + E^ 2 ]) 

[E(2/i + 3«:5 2 )] 2 

where 

Fi = fi[a i+1 - 2ai + Oj-i] - n[35 2 ai + Q5iV 2 } . 
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The second derivative, T, introduces the Lagrange multiplier through the equations of mo- 
tion: 

{ Oj = Vi = Fi - AViT } . 

Though complicated, the resulting equation for A is linear, and gives an explicit expression 
for the Lagrange multiplier that stabilizes T: 

•> -^Newton 

= E V,T ([2/x(F i+1 - 2F { + FUj - 6d%F\] - [E ^ 2 E 6*^]) ' 

Because T cannot change, both T = Tq and T = need to be properly specified in the 
initial conditions. A convenient way to do this is to choose the velocities equal to zero so 
that f = 0. 

A typical profile using two configurational thermostats is shown in Figure 6. Just as in 
the Hoover-Leete isokinetic case there is no equilibration of the constrained temperature 
(here configurational). The unconstrained kinetic temperature equilibrates throughout the 
system to a value dependent upon the initial conditions. Evidently the Landau-Lifshitz 
constrained configurational thermostat is quite useless for nonequilibrium work. 

X. CONCLUSIONS 

The results obtained with the Gaussian isokinetic and Nose-Hoover thermostats are quite 



QQflQ. 



For these 



consistent with the past history of their use over the last quarter century 
familiar thermostats there was no problem in reaching a nonequilibrium steady state with 
temperature profiles and heat fluxes very similar to those expected from macroscopic hydro- 
dynamics using Fourier's law. On the other hand, the results obtained with purely Hamilto- 
nian mechanics, using either the Hoover-Leete constrained isokinetic or the Landau-Lifshitz 
constrained configurational thermostats were both unexpected and thought-provoking. Both 
Hamiltonian approaches stably constrained their target temperatures in the two reservoir 
regions. Nevertheless temperature differences led to no net heat flow and the computed 
temperature profiles bore no resemblance to the predictions of hydrodynamics. The kinetic 
and configurational temperatures, though unequal in the constrained reservoirs, 

Tkin 7^ T con [reservoirs] , 
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were equal (to yet a third nearly-constant temperature) in the unconstrained bulk Newtonian 
regions: 

T kin = T con [bulk] . 

These results emphasize the fact that nonHamiltonian dynamics is absolutely necessary to 
the realistic modeling of nonequilibrium systems. A direct proof /understanding of the failure 
of the Hamiltonian systems to show a nonvanishing heat flux would be most welcome. 
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Figure Captions 

1. Temperature profile according to Fourier's law, with mean reservoir temperatures of 
0.26 (-200 < x < -100) and 0.24 (0 < x < 100). 

2. Temperature profile for a conventional two-temperature Nose-Hoover simulation, as de- 
scribed in the text. 

3. Temperature profile for a conventional two-temperature Gaussian isokinetic simulation, 
as described in the text. 

4. Summed-up values of the heat transferred to the hot and cold Gaussian isokinetic reser- 
voirs for the simulation of Figure 3. 

5. Temperature profiles (kinetic and configurational) for a simulation using constrained 
Hoover-Leete isokinetic reservoirs. The timestep is dt = 0.005, with a total elapsed time of 
10 6 divided equally between equilibration and averaging portions. 

6. Temperature profiles (kinetic and configurational) for a simulation using constrained 
Landau-Lifshitz isoconfigurational reservoirs. The timestep is dt = 0.002, with a total 
elapsed time of 400, 000 divided equally between equilibration and averaging portions. 



14 



[1] Wm. G. Hoover, Computational Statistical Mechanics (Elsevier, New York, 1991). 

[2] Wm. G. Hoover, Computer Simulation, Time Reversiblity, and Chaos (World Scientific Pub- 
lishing, Singapore, 1999 and 2001). 

[3] W. T. Ashurst and W. G. Hoover, "Argon Shear Viscosity via a Lennard-Jones Potential with 
Equilibrium and Nonequilibrium Molecular Dynamics", Physical Review Letters 31, 206-208 
(1973). 

[4] C. P. Dettmann and G. P. Morriss, "Hamiltonian Formulation of the Gaussian Isokinetic 

Thermostat", Physical Review E 54, 2495-2500 (1996). 
[5] D. J. Evans, W. G. Hoover, B. H. Failor, and B. Moran. "Nonequilibrium Molecular Dynamics 

via Gauss' Principle of Least Constraint", Physical Review A 28, 1016-1021 (1983). 
[6] W. G. Hoover, "Atomistic Nonequilibrium Computer Simulations", Physica 118A, 111-122 

(1983). 

[7] W. G. Hoover, "Mecanique de Nonequilibre a la Californienne" , Physica 240A, 1-11 (1997). 
[8] T. M. Leete, "The Hamiltonian Dynamics of Constrained Lagrangian Systems" [M S Thesis, 

West Virginia University, 1979]. 
[9] S. Nose, "Constant Temperature Molecular Dynamics Methods", Progress of Theoretical 
Physics Supplement 103 (Molecular Dynamics Simulations, S. Nose, Editor), 1-46 (1991). 
[10] L. V. Woodcock, "Isothermal Molecular Dynamics Calculations for Liquid Salts", Chemical 

Physics Letters 10, 257-261 (1971). 
[11] W. T. Ashurst, "Dense Fluid Shear Viscosity and Thermal Conductivity via Nonequilibrium 
Molecular Dynamics" [Ph D Dissertation, University of California at Davis/Livermore, 1974]). 
[12] L. D. Landau and E. M. Lifshitz, "Statistical Physics", equation 33.14 (McGraw-Hill, New 
Jersey, 1958). 

[13] J. G. Powles, G. Rickayzen, and D. M. Heyes, "Temperatures: Old, New and Middle- Aged" , 

Molecular Physics 103, 1361-1373 (2005). 
[14] K. Aoki and D. Kusnezov, "Nonequilibrium Steady States and Transport in the Classical 

Lattice 4 Theory", Physics Letters B 477, 348-354 (2000). Conductivity is 2.72T- 1 - 38 . 
[15] K. Aoki and D. Kusnezov, "Violations of Local Equilibrium and Linear Response in Classical 

Lattice Systems", Physics Letters A 309, 377-381 (2003). Conductivity is 2.8T" 1 - 32 . 



15 



Continuum Profile 



0.27 



0.25- 



0.23 




-200 







+200 



Gauss Temperature Profiles 



0.36 








0.32- 



0.30 



-200 







+200 



Gauss Heat Flux 




< Time < 200,000 



HL Temperature Profiles 



0.55 



1 1 1 

0.50 -tA^t+J^p** ^ > ^»swj ii, ^ w »' 

0.45- 
0.40- 
0.35- 
0.30- 
q 25 



0.20 



-200 



+200 



LL 250K Profiles 



0.35 



HOT(x) 



0.25- 



COLD(x) 




0.15 



Landau-Lifshitz T(x) 

i i i 



-200 







+200 



